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One of the most important effects of strong electron correlation in high-r c cuprates is the 
formation of Zhang-Rice singlets. By fully accounting for the quantum correlation effect of 
Zhang-Rice singlet formation, we show that a topological spin texture, skyrmion, is created 
around a Zhang-Rice singlet in the single-hole-doped Cu0 2 plane. The skyrmion picture 
provides a natural connection between the antiferromagnetic correlation and the doping con- 
^ 1 centration x. 
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1. Introduction 

It has been over two decades since the discovery of hole-doped high-critical-temperature 
(high-T,) superconductivity in copper oxides. However, despite the many novel pictures 

"\t" ■ 

presented, the mechanism behind the phenomenon remains elusive. One of the difficulties 
in explaining high-T^ cuprates is the presence of a strong electron correlation, which is re- 
sponsible for the Mott insulating state 2 in the parent undoped compound. Understanding this 
strong correlation in the doped Mott insulator would be the key to unveiling the mechanism 
behind the superconductivity. 

The Mott insulating state of the parent compound is a type of charge-transfer insula- 
tor. 3 The basic electronic structure of the Cu0 2 plane, which is the key structure in high-r c 
cuprates, is described by the d-p model. 4 An important observation was made by Zhang and 
Rice 5 about the description of a doped hole. Owing to the presence of a Cu-0 hybridization 
term, there is a strong antiferromagnetic (AF) interaction J K between an O p-hole and a Cu 
J-hole. When one considers combinations of four O-hole states around a Cu site, the sym- 
metric state with respect to the Cu site interacts with the Cu spin through an enhanced AF 
interaction. This strong correlation leads to a singlet state — the so-called Zhang-Rice singlet. 
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By restricting the Hilbert space where every doped hole forms a Zhang-Rice singlet and by 
excluding doubly occupied sites, the d-p model is reduced to the single-band t-J model, 
which has been studied extensively. 6 

The original analysis by Zhang and Rice was based on a small-cluster Hamiltonian con- 
sisting of a single Cu site and four surrounding O sites. However, if one takes a different point 
of view, the problem of describing Zhang-Rice singlet formation is equivalent to the prob- 
lem of describing a singlet state in a Mott insulator. However, here, in contrast to that in the 
cluster analysis, many-body effects should be considered. This kind of problem is generally 
quite difficult. In fact, describing a singlet state in a metal has been one of the central prob- 
lems in condensed-matter physics since its discovery. 7 8 In particular, a doped hole does not 
just form a Zhang-Rice singlet state, but has a certain amplitude of hopping to neighboring 
sites. Another nontrivial effect is that of the molecular fields created by nearest-neighbor Cu 
spins. In the presence of molecular fields, a singlet state changes from the pure singlet state 
to a state having a finite induced moment. In this study, we show that both of these effects 
play an important role in forming a topological spin texture, skyrmion, around a Zhang-Rice 
singlet. The physical picture that a skyrmion is formed around a Zhang-Rice singlet was pro- 
posed by the present author in ref. 9, and the mechanism of d^.^-wave superconductivity 
based on skyrmions was proposed in ref. 10. However, the microscopic mechanism of form- 
ing a skyrmion has not been established yet. Here, we demonstrate that a skyrmion is formed 
around a Zhang-Rice singlet on the basis of a microscopic model for high-r c cuprates. 

The skyrmion spin texture is a solution of the 0(3) nonlinear cr model (NLcrM). 11 The 
low-energy effective theory of the two-dimensional AF Heisenberg model is described by 
NLcrM: 12 



1 NLcrM 



- fdx fdy f d(c m f) \- (d t n) 2 - (d x ii) 2 - (d y vtf 



(1) 



Here, the continuous unit- vector field n = n(r, t) in a two-dimensional space r = (x, y) 
represents the local direction of the Neel order parameter, c sw is the spin-wave velocity, and g 
is the coupling constant. If we consider a static configuration of n, the states can be classified 
into homotopy sectors 13 characterized by the topological charge 

Q = J dx J dyp tc (x,y), (2) 

with the topological charge density 

Prc(r) = • (d x ii x d y n - d y n x d x n). (3) 

It is important to note that 13 any states in a given sector cannot move into another sector by 
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continuous deformation. The lowest energy in each sector is given by 



11 



(4) 



The skyrmion is the lowest-energy state in the Q = 1 sector. The explicit form is 

/ 2Ax 2Ay r 2 - A 2 \ 




(5) 



with r = sjx 2 + y 2 being the radial coordinate in space. 



In the early days of high-r c research, a skyrmion spin texture was conjectured 14 to corre- 
spond to the spinon excitation in the resonating- valence-bond theory. 15 When a Hopf term is 
added to S nlo-m, skyrmions can have fermionic statistics. 16 However, it was shown later that 
there is no Hopf term in the two-dimensional AF Heisenberg model. 1719 Thus, the possibility 
of identifying a skyrmion with a spinon was ruled out. 

Since a skyrmion is a metastable solution of NLcrM, there must be additional terms nec- 
essary to stabilize it. Such terms are inferred by scaling analysis where (x,y) —> (ax, ay). One 
can see that the energy of the NLcrM is scale-invariant, which is consistent with the skyrmion 
energy being independent of A. If there is a term that scales as 1 /a, it is possible to have a sta- 
ble skyrmion. For example, in the v = 1 quantum Hall system, one finds that the long-range 
Coulomb interaction scales as I /a. A skyrmion is stable provided that the Zeeman energy is 
sufficiently small. 20 In a doped antiferromagnet, Shraiman and Siggia pointed out 21 that the 
AF long-range order in the undoped system gives way to a non-collinear spiral state after hole 
doping. This spiral correlation is induced by a hole hopping term. The competition between 
the AF correlation and the spiral correlation leads to a modified hopping term that has an 
SU(2) matrix form in spin space. The components of the SU(2) matrix are determined by 
the spin configuration along the hopping direction. 22 In the continuum, one can see that this 
hopping term scales as I /a. However, it turns out that just having this term is not sufficient to 
stabilize a skyrmion. 

In this paper, we fully account for quantum effects in the formation of Zhang-Rice sin- 
glets in the single-hole-doped Cu0 2 plane to show that a skyrmion spin texture is formed 
around the singlet. The crucial point is that we need both the presence of a Zhang-Rice sin- 
glet and the hopping effect favoring a spiral state. 

This paper is organized as follows. In §2, we introduce the model that describes the single- 
hole-doped Cu0 2 plane. In §3, we present the numerical diagonalization result that suggests 
that a half-skyrmion spin texture is formed around a Zhang-Rice singlet. In §4, we present 
the self-consistent skyrmion configuration. In §5, the physical consequences of the skyrmion 
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picture are discussed. Section 6 is devoted to the conclusion. 
2. Model 

We consider a doped hole in the Cu0 2 plane in a skyrmion spin texture background. In the 
single-hole-doped Cu0 2 plane, the system gains energy from Cu-0 hybridization that leads 
to an AF interaction between O and Cu holes. 5 Zhang and Rice 5 considered the combinations 
of the four O-hole states around a Cu ion, and showed that the singlet state between a Cu 
hole and a symmetric O-hole state with respect to the central Cu ion has the largest binding 
energy. By projecting out the other O-hole states, the t-J model was derived. 5 Although 
we assume that the correlation of forming this singlet state is the most important electronic 
correlation, we do not restrict ourselves to the Hilbert space where the singlet is fully formed. 
The Hamiltonian is 

+ f Z ^r(pWp e ) + jJ]SrSj, (6) 

Here, the indices £ and £' refer to O sites, and i and j refer to Cu sites. The hole hopping 
integral t is restricted to nearest-neighbor sites, and the components of cr = (cr x , cr y , cr z ) are 
the Pauli matrices. The exchange interaction J between Cu spins S ; is restricted to nearest- 
neighbor sites. The operator p\ a creates an O (2p x , 2p y ) hole with spin cr at site i. In order to 
describe Zhang-Rice singlet formation, the Cu spin at the center is described by the quantum 
spin 

So = i(J W ), (7) 

with d\ = (Jq T , d^). The operator creates a Cu 3d x 2_ y 2 hole with spin cr at the central Cu 
site. The second term of eq. (6) describes the interaction of forming a Zhang-Rice singlet. 
The operator if/ consists of four O-hole states around the quantum Cu spin at the center and 
is given by 

MO) 

If we take a strong J K limit, the Zhang-Rice singlet is fully formed and eq. (6) is reduced to 
the t-J model. However, here, we do not take this limit. 

We assume a skyrmion configuration for S 7 except at the origin (j = 0). By representing 
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the coordinate vector of the Cu site j by R y , S 7 is given by 

Sj = (-iySn sk (Rj,A), (9) 

where n sk is defined by eq. (5). For S , we set S = 1 /2. (If we take into account the effect of 
quantum fluctuations, S is somewhat reduced. However, we neglect this for simplicity.) We 
use NLcrM in evaluating the energy associated with the last term in eq. (6). The energy of a 
spin texture in NLcrM is given by eq. (4). The salient feature of this energy formula is that 
the energy is independent of the skyrmion size A. The resulting spin texture can be different 
from the originally assumed skyrmion configuration because of the presence of a quantum 
Cu spin at the origin. However, the static Cu spin configuration is still characterized by the 
topological charge Q, and the energy is given by eq. (4) as long as we consider the lowest 
energy state in the homotopy sector with Q. Since eq. (6) is independent of the skyrmion size 
A, we may focus on the energy of the hole and the quantum spin S and their interaction with 
other Cu spins. Therefore, the Hamiltonian eq. (6) is reduced to the following form: 

*H = ~t J] (PvtrPto- + h -c) + Jk (dlo-do) ■ (^JotAo) 

+y Yj s r(pl (T Pe) + i^]( d l ' d o)-^j + E NL(T M- (10) 

K*o),t <oj> 

In the second last term, we consider the interaction between S and its nearest neighbor Cu 
spins. In the analysis below, we denote energies in units of t. 

3. Skyrmion Creation around Zhang-Rice Singlet 

The Hamiltonian eq. (10) is exactly diagonalized in the Hilbert space where each of the 
central Cu-hole and the O-hole states is singly occupied. For J, we assume J/t =1/3, which 
is a standard value taken from the literature. 6 The A dependence of the ground-state energy 
is shown in Fig. 1. The ground-state energy exhibits a clear energy minimum for J K > 1. In 
particular, the minimum is located aX A - a, and this value is insensitive to an increase in J K 
for J K > 1.5. 

Now, we estimate J K . From the second-order perturbation theory with respect to the Cu-0 
hybridization term, 5 J K is given by 

Jk = -^t + 7^' (ID 
U d - A U p + A 

where td p is the wave-function overlap of Cu and O holes. U p is the on-site Coulomb repulsion 
at an O site and Ud is the on-site Coulomb repulsion at a Cu site. The parameter A is given 
by A = s p - e d , with e p and e d being the atomic energies of the O and Cu holes, respectively. 
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Fig. 1. (Color online) Skyrmion core size A dependences of the ground-state energy of the doped hole Hamil- 
tonian eq. (6) for different values of Jk- The energy without a Zhang-Rice singlet for Jk — 2 is also shown. The 
energies are measured from their values at A — 0. In describing the CuC<2 plane, we employed a 21 x 21 square 
lattice for Cu sites and a 12 x 12 square lattice for O sites with open-boundary conditions. 



Using the parameter set taken from constrained local density approximation calculations, 
we find J K 2. Therefore, a skyrmion with a core size of A - a is formed around the Zhang- 
Rice singlet. Since the core state is the Zhang-Rice singlet state, where there is no topological 
charge density, this result suggests that a half- skyrmion with \Q\ = 1/2 is created around the 
Zhang-Rice singlet. It is important to note that the presence of the Zhang-Rice singlet is 
crucial to stabilization of the skyrmion. Thus, if one replaces the quantum Cu spin with the 
classical spin, the minimum disappears, as shown in Fig. 1. 

The location of the energy minimum at A - a can be explained by a simple calculation. 
Let us consider a cluster consisting of a single Cu site at (0, 0) and four surrounding O sites at 
(+a/2, 0), (0, ±a/2). We assume that there are molecular fields at the O sites created by other 
Cu spins at (±a, 0) and (0, +a). Then, the exact diagonalization of this cluster Hamiltonian 
is carried out analytically. (The details of the calculation is presented in Appendix .) For the 
radial configuration of the molecular fields, which is equivalent to the skyrmion configuration 
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with A = a, we find that the energy is given by 

3J K + It 1 
2 



E r = -- 



-^(3J K + 2t) 2 + S 2 J 2 K . 



Meanwhile, for the AF configuration of the molecular fields, we find 



= -J K -2t- ^27-^) 2 52 + 47 |_ 



(12) 



(13) 



Using the parameters for the Cu02 plane, we find E r < E AF as shown in Fig. A-2. This can be 
simply understood as follows: In order to maximize the energy gain by creating the Zhang- 
Rice singlet state, the system prefers to have the molecular fields cancel each other. This is 
achieved for the radial configuration. By contrast, the molecular field effect is maximized for 
the AF configuration. 

The instability of the homotopy sector with Q = 0, which corresponds to the AF state, is 
understood in the following way as well. In order to clarify the physical picture, we rewrite 
the NLcrM by introducing complex fields z a (r, through 

n(r,0= Yj 4(r,0ov,^(r,0, (14) 

<r,o-'=U 

and we obtain the CP 1 model 13 



S C p l = ^ j dx j ' dy J d(c SVi t) ^ 



— d t - ia Q I Za (r, t) 



- \{d x + ia x ) Z(T (r, Ol 2 - \(d y + ia y ) z a (r, t) 



with 



a = i V 4 (r, t) — - Za (r, t) , 
f-J d(c sw t) 

a= l -Y J K ( r > Vz<r (r, t) - [W a (r, t)] Za (r, t)} . 

<r=tl 

The equation of motion for the steady state is given by 

(V + iafza = 0. 

In the presence of the Zhang-Rice singlet, the boundary condition at the origin is 

Za(0) = 0. 



(15) 

(16) 
(17) 

(18) 
(19) 



In the AF configuration, there is no gauge flux, V x a = 0. Therefore, the equation of motion 
is reduced to 



V 2 Za = 0. 



(20) 
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Using the polar coordinate r and with x = r cos <p and y = r sin (p, let z a = f a (r) exp(i£(p). 
Here, I is an integer. The equation for f^r) is given by 

^ + I^_^ = . (21) 
dr 2 r dr r 2 

To be consistent with the assumption that there is no gauge flux, we must take t = 0. The 
equation is readily solved and we find 

/^oclnr. (22) 

Clearly, this solution is incompatible with the boundary condition eq. (19). 

In the case that there is a nonvanishing gauge flux, it is difficult to solve the equation of 
motion because it is a nonlinear equation with respect to z<r(r, t). However, it is possible to 
find the behavior of the solution near r = 0. If we assume z a = fA r ) exp(z'^>) with I ± 0, we 
find 

a = -~Y J fl^ (23) 

(T 

where = (- sin <f>, cos cf)). From eq. (18), after some algebra, we obtain 

iH£-?Z(i-*JW ™ 

cr' 

For r —> 0, the solution of this equation behaves as 

U ^ Cj { \ (25) 

with Co- as a constant. Note that this form is consistent with the boundary condition eq. (19). 
Since I ± 0, there is a gauge flux. The situation is quite similar to the problem of the point 
singularity in a two-dimensional superfluid. 24 The vortex in Zo-(r) turns out to be a skyrmion. 
In fact the topological charge density eq. (3) is expressed as 

Ptc(r) = (d x a y - d y a x ) , (26) 

in terms of the gauge flux. 13 Therefore, a nonvanishing gauge flux is a skyrmion characterized 
by the topological charge eq. (2). 

4. Self- Consistent Skyrmion Configuration 

In the previous section, we have presented the numerical result that suggests that a half- 
skyrmion spin texture is formed around the Zhang-Rice singlet. In this calculation, Cu spins 
except that at the origin are assumed to be a skyrmion spin texture given by eq. (9). In this 
section, we attempt to obtain the self-consistent spin configuration by fixing the directions of 
the four Cu spins surrounding the Zhang-Rice singlet. 
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The self-consistent Cu spin configuration is obtained by solving the mean field equations. 
At each Cu site, there are molecular fields created by the nearest-neighbor Cu sites and O- 
hole. In the self-consistent configuration, the Cu spin at each site is antiparallel to this molec- 
ular field. To describe the Cu0 2 plane, we employ a 15 x 15 square lattice for Cu sites and a 
10 x 10 square lattice for O sites under the open-boundary condition. As explained below, the 
doped hole is strongly localized around the Zhang-Rice singlet. Therefore, a large system 
is not necessary in describing doped hole hopping. According to the energy calculation in 
Fig. 1 and the cluster analysis above, we fix the four Cu spins around the center. As shown 
in the self-consistent spin configurations in Fig. 2, the difference between the skyrmion and 
the anti-skyrmion is only that the ^-components of the spins are reversed. Furthermore, as 
shown in Fig. 2(c), the doped hole is strongly localized around the quantum Cu spin. In the 
absence of the Zhang-Rice singlet, the doped hole is not localized but extended. Therefore, 
the formation of a skyrmion is not sufficient to localize the doped hole. 

The physically important quantity is the topological charge density shown in Fig. 2(d). 
Note that rotation of all spins about the z-axis through the same angle does not change the 
topological character. However, the topological charge changes sign for the anti-skyrmion 
case. It is important to note that the topological charge density is considerably extended when 
compared with the doped hole wave function. Since the topological charge distribution pro- 
duces a magnetic-field-like effect through the Berry phase, 25 the non-vanishing topological 
charge density carried by the doped hole produces intriguing effects, as will be discussed be- 
low. The sum of the topological charge is 0.320, whereas the expected value is one-half. This 
discrepancy is probably associated with the approximation in describing Cu spins by classi- 
cal spins. In fact, the calculation above overestimates the AF correlation that reduces the total 
topological charge. Describing Cu spins using more suitable fields is a future problem. 

The existence of a stable skyrmion spin texture in a Li-doped system was suggested by 
Haas et al. 26 in the t - J model with one hole on a 13-site cluster. In the cluster, the central 
Cu site is replaced by Li. The exact diagonalization of the cluster system under the boundary 
condition imposed by the skyrmion configuration shows that the total energy has a minimum 
at A - 2a. In their calculation, the doped hole cannot reside at the center site and is bounded 
by the attractive Coulomb potential of Li. Although the system is not equivalent to our system, 
the mechanism of stabilizing a skyrmion is quite similar. In the model in ref. 26, the doped 
hole is localized owing to the attractive Coulomb potential. In our system, however, the doped 
hole is localized by the formation of the Zhang-Rice singlet. The stability of a skyrmion is 
also demonstrated in a 16-site cluster 27 with hole motion being restricted to four sites around 
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Fig. 2. (Color online) Self-consistent Cu spin configuration in the case of a skyrmion (a) and an anti-skyrmion 
(b). The dashed arrows represent the doped hole spin, which is scaled by a factor of 10. (c) Doped hole density 
and (d) topological charge density for the skyrmion (a). 



the center. 

5. Discussion 

Now, we explore the physical consequences of the skyrmion spin texture created by a 
doped hole. The most direct consequence is that the magnetic properties are straightforwardly 
connected with the hole doping concentration x. As revealed by neutron scattering measure- 
ments, 28 the inverse of the AF correlation length % AF is well described by the simple formula 

k(x, T) = k(x, 0) + k(0, T), (27) 

where T is the temperature. The inverse of the temperature-independent term k(x, 0) is given 
by 

1 a 

-—^r = -p- (28) 
k(x, 0) 
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The right-hand side of this equation is the average separation between the doped holes. Mean- 
while, k(0, T) is described by the formula for NLcrM. 29 The latter is associated with the under- 
lying AF correlation. When each doped hole creates a skyrmion, local spin disorder is created 
by the skyrmion. Therefore, the average hole separation leads to a temperature-independent 
component of %af- Another important observation is that the skyrmion character is lost at 
approximately around x ~ 0.3, where the average hole separation is ~ 2a. In the vicinity of 
this critical doping level, each doped hole is unable to create a skyrmion with A = a and so 
all skyrmion features are lost for x > 0.3. 

As we have seen above, the doped hole is localized by the formation of the Zhang-Rice 
singlet. Therefore, Nagaoka ferromagnetism 30 does not occur and the system is insulating in 
a slightly doped regime. The energy band of doped holes forms when the overlap of adjacent 
doped hole wave functions is not negligible. The bandwidth W is proportional to the wave 
function overlap. 

As mentioned above, the presence of a skyrmion leads to a Berry phase effect. The 
skyrmion produces a magnetic-field-like effect for doped holes. When a doped hole enters 
a region of skyrmion spin texture created by another doped hole, then the Lorentz force effec- 
tively acts between the two holes. This interaction leads to a chiral pairing state. 10 Since there 
are skyrmions and anti-skyrmions, two chiral pairing states coexist with opposite chiralities. 
The analysis of the gap equation suggests that this type of interaction leads to d x 2_ y 2 -wave 
superconductivity. 1031 

In addition, skyrmion dynamics has a nontrivial contribution to hole dynamics. The 
skyrmion itself has its own dynamics, which is not identical to bare hole dynamics. 9 The 
doped hole moves with a constraint of keeping the skyrmion as long as W < J. The un- 
derlying Fermi surface, which appears when the skyrmion is suppressed, is partially gapped 
because the skyrmion energy does not necessarily vanish along it, which results in a trun- 
cated Fermi surface. 31 For this pseudogap phenomenon, the characteristic temperature is the 
binding energy of the skyrmion, which is ~ 0.04? from Fig. 1. The order of magnitude for 
this value is the same as that of the pseudogap temperature evaluated in materials with high 
T c cuprates. 32 However, for elaborate investigations of the pseudogap, a better description is 
necessary. In this regard, it is worth pointing out that the Hamiltonian describing the dynam- 
ics of the composite object comprising the doped hole and skyrmion is similar to that for 
quasiparticles in the J-density wave state. 33 However, the physical picture is quite different. 

Here, we comment on electron-doped high-T c cuprates. The superconducting state of 
electron-doped cuprates has a rf-wave pairing symmetry, 34 which is the same as that in the 
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hole-doped case. From this observation one might expect that the mechanisms of supercon- 
ductivities are the same for electron-doped and hole-doped cuprates. However, the super- 
conducting transition temperature of electron-doped cuprates is quite low compared with 
that of hole-doped cuprates with high transition temperatures. Furthermore, there are quali- 
tative differences between electron-doped and hole-doped cuprates. Among others, the rela- 
tionships between superconductivity and antiferromagnetism are strikingly different. In hole- 
doped cuprates, only 2-3% hole doping is sufficient to kill AF order, while we need about 
20% electron doping to suppress AF order in electron-doped cuprates. In the phase diagram 
in the plane of doping concentration and temperature, the superconducting phase is clearly 
separated from the AF phase in hole-doped cuprates, while the superconducting phase and 
the AF phase sit side by side in the electron-doped cuprates. 35 Therefore, it is quite natural 
to surmise a different mechanism for superconductivity works in electron-doped cuprates. If 
the mechanism is different, it is quite hard to describe such a difference in the single-band 
Hubbard model. Meanwhile, the difference is simply associated with the presence or absence 
of skyrmions in the skyrmion picture. In electron-doped cuprates, Zhang-Rice singlets are 
irrelevant because electron carriers are introduced at Cu sites and there is no singlet corre- 
lation between electron carriers and Cu J-electrons. Therefore, a skyrmion is not formed in 
electron-doped cuprates. 

Finally, we propose a "smoking-gun" experiment to establish the presence of skyrmions 
in high-r c cuprates. In the doping range of 0.02 < x < 0.055, incommensurate inelastic 
magnetic peaks are observed in neutron- scattering experiments. 36 37 In the skyrmion picture, 
these peaks are associated with the AF skyrmion lattice. 31 Recently, real-space observations 
of a two-dimensional skyrmion lattice have been reported in Feo.5Coo.5Si. 38 This experimental 
technique can be applied to finding the skyrmion lattice in high-r c cuprates as well. 

6. Conclusions 

To conclude, we have shown that a skyrmion spin texture is created around a Zhang-Rice 
singlet in the single-hole-doped Cu0 2 plane. The evolution of the AF correlation upon hole 
doping is naturally described by the skyrmion picture. A doped hole with a skyrmion creates 
the Berry phase effect, which can lead to J-wave superconductivity. The skyrmion picture 
also seems to be consistent with some other aspects of high-r c cuprates. 
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Fig. A'l. (Color online) Molecular fields Sj (j = 1,2,3,4) at Cu sites for the (a) AF configuration and (b) 
radial configuration. The circle at the center denotes the quantum Cu spin. The O sites are labeled as 1,2, 3, and 
4. 
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Appendix: Derivation of Eqs. (12) and (13) 

In this Appendix, we present the derivation of eqs. (12) and (13). We consider a cluster 
that consists of one quantum Cu spin sitting at the center, four O sites, and four classical Cu 
spins Sj, as shown in Fig. A-l. For convenience, we introduce the following field: 

y^ = \ exp ^ iq ~ ^ p<a ~' ( A *^ 

£=1,2,3,4 

where q = 0, n/2, n, 3n/2. In terms of these fields, the Hamiltonian of the cluster is given by 
^ = 2t J] Hkr&r + rLrna) + Jk (4<rd ) ■ (rJoTo) 

O" 

;'. p, i j 
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where (Mj) pq = (lP) pj U jq with 



U -2 



1 1 



-1 

-i 



-i 



1 1 
1 -1 



•1 1 1 
i 1 -1 



(A-3) 



For the AF configuration of S ; (Fig. A- 1(a)), we take Sj = Se z , where e z is the unit vector 
in the z-direction. The third term in eq. (A-2) has the following simple form: 



S Jr ^ 



y\o- z y p . 



(A-4) 



In the Hilbert space, where each of the central Cu states and O states is singly occupied, the 
Hamiltonian is diagonalized analytically, and we find that the ground-state energy is given by 
eq. (13). 

For the radial configuration of Sj [Fig. A-l(b)], we take Si = Se x , S2 = Se y , S 3 = -Se x , 
and S 4 = -Se y . Here, e x and e y are the unit vectors in the x- and v-directions, respectively. In 
this case, the third term in eq. (A-2) has the form 
SJ K 



(A-5) 



(A-6) 



where cr ± = (cr x + icr y )/2. Noting that this term does not contain y ncr or y^, we find that the 
Hilbert space is divided into subspaces. Explicitly considering spin states of the Cu spin state 
and the O hole state, the problem is reduced to diagonalize the following 4x4 matrix: 

-2t-J K 2J K s -f- 
2J K -2t-J K s -f 
^ 

^00 
By taking the lowest eigenenergy of this reduced Hamiltonian, we obtain eq. (12). 

The energies of eqs. (12) and (13) are shown in Fig. A-2 as functions of 5. Here, we 
assume J K /t = 2 and J/t = 1 /3. In the entire range of S , the energy of eq. (12) is lower than 
the energy of eq. (13). We find that E r < E AF for J K /t > 0.66. For J K /t < 0.66, E r > E AF in 
the entire range of S . 
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Fig. A-2. (Color online) E r and E AF as functions of the spin S for J K /t - 2. 
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